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ABSTRACT 

The ecologically important cyanobacterium 
Prochlorococcus possesses the smallest genome 
among oxyphototrophs, with a reduced suite of 
protein regulators and a disproportionately high 
number of regulatory RNAs. Many of these are 
asRNAs, raising the question whether they modulate 
gene expression through the protection of mRNA 
from RNase E degradation. To address this ques- 
tion, we produced recombinant RNase E from 
Prochlorococcus sp. MED4, which functions opti- 
mally at 12 mM Mg 2+ , pH 9 and 35 C. RNase E 
cleavage assays were performed with this recom- 
binant protein to assess enzyme activity in the 
presence of single- or double-stranded RNA sub- 
strates. We found that extraordinarily long asRNAs 
of 3.5 and 7 kb protect a set of mRNAs from RNase E 
degradation that accumulate during phage infec- 
tion. These asRNA-mRNA duplex formations 
mask single-stranded recognition sites of RNase 
E, leading to increased stability of the mRNAs. 
Such interactions directly modulate RNA stability 
and provide an explanation for enhanced transcript 
abundance of certain mRNAs during phage infec- 
tion. Protection from RNase E-triggered RNA 
decay may constitute a hitherto unknown regulatory 
function of bacterial c/s-asRNAs, impacting gene 
expression. 

INTRODUCTION 

Endoribonuclease E (RNase E) initiates the decay of most, 
if not all, mRNAs in prokaryotes by endonucleolytic 
cleavage followed by polyadenylation and/or exonucleo- 
lytic degradation (1). Therefore, RNase E is a crucial regu- 
lator of the global rate of RNA decay. The half-life of 



RNA is an important factor in the regulation of gene 
expression, directly influencing the translation rates of 
their encoded proteins and hence impacting the rapidity 
of the response of an organism to changing environmental 
conditions. Several studies have shown that stress condi- 
tions, such as the presence of antibiotics, nutritional stress 
and transitions in growth phase, cause a dramatic change 
in the rate of RNA turnover for a subset of genes (2-4). 

The presumably best characterized RNase E homologue 
is that of Escherichia coli, which is organized in a multimer 
enzyme complex known as the degradosome (5). The 
degradosome consists of polynucleotide phosphorylase, 
RNA helicase B and the glycolytic enzyme enolase (5-8). 
Degradosome structures have not been identified in 
cyanobacteria or spinach chloroplasts, though all compo- 
nents necessary for the assembly of the degradosome 
are present in cyanobacteria, including Synechocystis 
PCC6803 and Prochlorococcus sp. (9-11). However, 
cyanobacterial RNase E (me) genes of the species 
Synechocystis PCC6803 and Prochlorococcus sp. lack the 
carboxy-terminal sequence necessary for the assembly of 
the degradosome in E. coli (12). Nevertheless, similar and 
important functions of RNase E in both E. coli and 
Synechocystis PCC6803 can be assumed, as the disruption 
of rne is lethal in both organisms (11). The role of 
degradosome structures is still under debate because the 
deletion of the carboxy terminus in E. coli (associated with 
a loss of the degradosome structure) does not impact cell 
viability (11). According to the lack of the C-terminal 
extension, a feature of RNase G (Caf A family), and 
sequence homology to RNase E, cyanobacterial homo- 
logues have been grouped into the RNase E/G family (13). 

Several recent studies have been published on whole- 
transcriptome expression analyses, revealing that a signifi- 
cant portion of the transcriptome is expressed on the anti- 
sense strand not only in Eukaryotes (14-17) but also in 
Bacteria and Archaea, where data indicate that between 
3% to more than half of all genes are associated with 
asRNAs (18-24). In light of RNase E enzyme 
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specificity — recognizing single-stranded A/U-rich se- 
quences (25,26) — one wonders if the potential formation of 
asRNA-mRNA duplexes affects RNase E activity. In this 
context, in a recent microarray study on whole-genome 
expression of the podovirus P-SSP7 and its cyanobacterial 
host Prochlorococcus MED4, we found that transcript 
abundance of 41 host genes, including RNase E, increased 
during lytic infection, while the vast majority of host 
mRNAs declined in abundance (27). We suggested that 
the induction of host RNase E transcript may help the 
phage by generating a source of nucleotides that (after 
nucleotide reduction by the phage-encoded ribonucleotide 
reductase) are utilized by the phage for its own replication 
(27). Those results raised the question of how these 41 host 
mRNAs are protected from degradation in light of the 
enhanced transcription of RNase E. 

Here, we set out to test the hypothesis that in the ecologic- 
ally important marine cyanobacterium Prochlorococcus, 
asRNA-mRNA duplex formation masks single-stranded 
recognition sites of RNase E, leading to increased stability 
of mRNAs and their asRNAs. Therefore, using purified 
recombinant Prochlorococcus RNase E from E. coli, we 
assessed RNase E activity in the presence of single- or 
double-stranded RNA substrates and showed that 
asRNA-mRNA duplexes protected mRNA from RNase E 
degradation, enhancing mRNA stability of transcripts 
accumulated during phage infection. 

MATERIALS AND METHODS 

Construction of vne expression plasmid pQE70PMM1501 

Oligonucleotides used for PMM1501 cloning were 
PMM1501_C_K12 fw and PMM1501_C_K12 rev. 
Lowercase letters indicate Pael or BamHI site 
(Supplementary Table SI). Both oligonucleotides were 
optimized for E. coli codon usage (Kazusa Codon Usage 
Database) without impairing the PMM1501 wild-type 
amino acid composition. The resulting coding sequence 
of PMM1501 was cloned in pQE70 (QIAGEN, Germany). 

Overexpression and purification of recombinant RNase E 
from Prochlorococcus MED4 

Escherichia coli TOP10 F' served as the host strain for 
construction and propagation of pQE70-PMM 1 50 1 . 
Expression of recombinant PMM1501 protein was 
carried out in E. coli ArcticExpress (DE3) RIL 
(Stratagene, USA). Strains were grown in Luria broth, 
supplemented with ampicillin (lOOug/ml), streptomycin 
(75 |ig/ml) and gentamicin (20 ug/ml), if required. An over- 
night culture was diluted in 2 1 of Luria broth and grown 
at 37°C to an OD 600 of 0.5-0.8. Subsequently, heter- 
ologous expression was induced with isopropyl 
(3-D-l-thiogalactopyranoside (IPTG) (0.1 mM final con- 
centration), and the culture was transferred to 10°C and 
grown for an additional 48 h. Thereafter, cells were har- 
vested (4°C, 6000 rpm, lOmin) and resuspended in buffer 
A [50 mM NaH 2 P0 4 , 300 mM NaCl, 20 mM imidazole, 
(pH adjusted to 8.0 using NaOH)], containing protease 
inhibitor PMSF (0.2 mM). Disruption of cells was 
carried out on ice by sonication. Cell debris was 



removed by centrifugation at 13 000 rpm for 30min at 
4°C, and the supernatant was collected. C-terminally 
His-tagged MED4 RNase E was purified under native 
conditions by Ni 2+ nitrilotriacetic acid affinity chromatog- 
raphy (QIAGEN, Germany) according to the manufac- 
turer's instructions (QIAGEN handbook, fifth edition). 
All steps were performed at 4°C, but the pH of all 
buffers was adjusted at room temperature. The Ni-NTA 
matrix size was adjusted to 500 ul. Following equilibra- 
tion with 10 volumes of buffer A, the column was 
charged with supernatant and subsequently washed with 
60 volumes of buffer B [50 mM NaH 2 P0 4 , 500 mM NaCl, 
50mM imidazole, 0.1% (v/v) Triton X-100, 8% (v/v) 
glycerol (pH adjusted to 8.0 using NaOH)]. Adsorbed 
peptides were eluted by ligand exchange employing a 
step gradient of 250 and 500 mM imidazole in buffer con- 
taining 50 mM NaH 2 P0 4 , and 300 mM NaCl, pH 8.0 
(adjusted with NaOH). Fractions containing specific 
endoribonuclease activity were identified by in vitro 
assays as described below (in vitro cleavage assay) and 
visualized by Coomassie-stained 10% SDS-PAA gels. 

Size exclusion chromatography 

Combined elution fractions of His-tag-purified MED4 
RNase E were size fractionated on a superdex 
200 10/300 GL column (GE Healthcare, Germany). The 
column was equilibrated with buffer A (see purification 
of recombinant protein) and calibrated using a gel filtra- 
tion standard, ranging from 1.35 to 670 kDa (Bio-Rad, 
Germany). Fractionation was analysed and documented 
with Unicorn 4.12 software (Amersham Pharmacia, 
Germany). 

In vitro synthesis of RNA 

In vitro transcription was carried out with T7 polymerase 
and PCR fragments containing a T7 recognition site as 
templates using MEGAshortscript™ High Yield 
Transcription KIT (Ambion, USA), according to the 
manufacturer's instructions. Residual DNA was 
removed with 2U TURBO DNase (Ambion, USA) for 
15min at 37° C. RNA was purified and concentrated by 
phenol-chloroform extraction and ethanol precipitation 
and finally quantified using a NanoDrop ND-1000 spec- 
trophotometer (PEQLAB biotechnology, Germany). Ten 
micrograms of purified in vitro transcript was subjected to 
electrophoretic separation on a 7M urea-6% polyacryl- 
amide gel. The RNA fragment, corresponding to full- 
length transcript, was excised and recovered from the gel 
slice by electro-elution for 45 min at 1 1 mA per glass tube 
(Model 422 Electro-Eluter; BIO-RAD, Germany). 
Enrichment of gel-purified RNA was conducted by subse- 
quent gel chromatography using RNA Clean & 
Concentrator™-25 according to the manufacturer's in- 
structions (Zymo Research, USA). If required, RNA 
was treated with tobacco acid pyrophosphatase or 
S'-polyphosphatase (Epicentre, USA) according to the 
manufacturer's instructions and subjected to an additional 
purification by phenol-chloroform extraction and ethanol 
precipitation (prior to gel purification). The following 
mRNAs and asRNAs were synthesized: PMM0684 
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RNA, PMM0685 RNA, PMED4_07431 RNA, PMED4_ 
07411 RNA, 3kb polycistronic mRNA and 3.5 kb 
asRNA. Oligonucleotides used for template generation 
are listed in Supplementary Table SI. 

In vitro cleavage assay 

For ribonuclease activity tests of duplex RNA, both 
mRNA and asRNA were diluted in H 2 0 in the same 
reaction vial, heated to 95°C for 3-5 min to break intrinsic 
secondary and tertiary structures, and then briefly chilled 
on ice. The reaction mixture was complemented with 
25 mM Tris-HCl (pH 8.0), 60 mM KC1, 5mM MgCl 2 , 
lOOmM NH 4 C1 and 0.1 mM DTT. Hybridization of com- 
plementary sequence regions was carried out at room 
temperature for 20 min. 

RNase E. In vitro cleavage was carried out for 15 min at 
30° C in a 5-ul reaction mixture containing 25 mM Tris- 
HCl (pH 8.0), 60 mM KC1, 5mM MgCl 2 , lOOmM NH 4 C1, 
0.1 mM DTT, 5% (w/v) glycerol, 0.2-1.3 pmol ssRNA or 
0.4-2.6 pmol dsRNA, and 7 pmol affinity-purified RNase 
E. Endoribonuclease activity was abolished by addition of 
1 |il 0.5 M EDTA and 1 volume loading buffer or by direct 
addition of 1 volume loading buffer. Finally, samples were 
heated to 95°C for 3-5 min and briefly chilled on ice before 
the cleavage reactions were analysed by employing elec- 
trophoretic separation on 7M urea-6% polyacrylamide 
gels for 2 h at 6 mA. 

RNase III and RNase Tl. Assay conditions were identical 
to those of the RNase E approach, with the following 
modifications. RNA digestion was carried out at 37°C 
for 15 min using 0.1 U RNase III (Epicentre, USA; from 
E. coli) and 0.5 U RNase Tl (Sigma, Germany; from 
Aspergillus oryzae), respectively, or according to experi- 
mental requirements. 

Fluorogenic cleavage assay 

A fluorogenic RNA oligonucleotide with a FAM tag and 
a BHQ-1 quenching tag containing the RNase E recogni- 
tion site of Synechocystis PCC6803 in the psbA gene 
leader sequence, psbA2-5 / UTR [5 / -Phos-A-A(2 / -0- 
Methyl)-C(2 / -0-Methyl)-A(2 / -0-Methyl)-U(FAM)-C(2 / - 
0-Methyl)-G(2 / -0-Methyl)-A(2 / -0-Methyl)-A-A-A-U- 
A-C-A-U-A-C(BHQ-l)-OH], was synthesized and 
PAGE purified by Purimex (Grebenstein, Germany). 
Digestion of the fluorogenic oligonucleotide results in a 
marked increase in fluorescence that can be monitored 
in real time under defined reaction conditions. Cleavage 
assays were carried out at 30°C in a 50-ul reaction mixture 
containing 25 mM Tris-HCl (pH 8.0), 60 mM KC1, 5mM 
MgCl 2 , 100 mM NH4CI, 0.1 mM DTT, 5% (w/v) glycerol, 
1 |il affinity-purified RNase E (not quantified) and 0.1 (iM 
(in combination with MED4 RNase E) or 0.06 uM (in 
combination with PCC6803 RNase E) fluorogenic RNA, 
or adjusted according to experimental requirements. To 
avoid evaporation during the course of procedure, 
samples were covered with 20 ul RNase-free immersion 
oil. Each cleavage experiment was assayed in technical 
triplicates. The cleavage reaction was monitored 



continuously by fluorometry using a Victor™ X 3 
multilabel plate reader (PerkinElmer; excitation at 
480 nm, emission at 520 nm). 

Growth conditions of Prochlorococcus MED4 

Cells were grown at 22° C in AMP1 or Pro99 medium (28) 
under 30 umol quanta m -2 s _1 continuous white cool light. 
For phage infection experiments, the podovirus P-SSP7 
was added at a ratio of three infective phage per cell 
(MOI = 3) such that the cell and phage concentrations 
were 10 8 cells/ml and 3 x 10 8 phage/ml, respectively. 
Prior to phage addition, the cells were growing exponen- 
tially with a generation time of 1.5 days as were the unin- 
fected control cells. 

RNA isolation 

Total RNA was extracted from cultures using the hot 
phenol method (29,30). Briefly, cells were harvested by 
centrifugation at 12 000 rpm for 10 min at 20°C. The 
pellet was resuspended in RNA resuspension buffer 
[10 mM sodium acetate (pH 5.2), 200 mM sucrose, 5mM 
EDTA]. Following the addition of half a volume of 
phenol, the suspension was incubated in a 60° C water 
bath for 15 min and subsequently incubated at room tem- 
perature for another 30 min. Half a volume of chloroform/ 
isoamyl alcohol (24:1) was added, and the homogenate 
was thoroughly mixed, incubated for another 1-10 min 
at room temperature and centrifuged at 6000 rpm 
for 5 min at 20°C. The upper aqueous phase was extracted 
once more with one volume phenol/chloroform/ 
isoamyl alcohol (25:24:1) and centrifuged for 5 min at 
6000 rpm (20° C). The aqueous phase was precipitated 
with 1/10 volume 3M sodium acetate, 3 volumes 100% 
ethanol and 0.2 ul glycogen at — 20°C overnight. 
Subsequently, RNA was concentrated by centrifugation 
at 13 000 rpm for 30 min (4°C) and resuspended in 
RNase-free H 2 0. 

In vitro synthesis and labelling of RNA 

Oligonucleotide probes. Ten picomoles of oligodeoxy- 
nucleotide were radiolabeled at the 5'-end with 10U T4 
polynucleotide kinase (Fermentas, Germany), 3 nmol 
[y- 32 P]ATP (Hartmann Analytic, Germany), 5U 
RiboLock RNase inhibitor and 1 x reaction buffer A 
(Fermentas, Germany) in a total volume of 20 ul. 
Following incubation at 37°C for 30 min, the reaction 
was stopped by the addition of 25 mM EDTA and de- 
naturation at 95°C for 5 min. The reaction mixture was 
diluted in H 2 0 to a total volume of 40 ul and subjected to 
separation on pre-packed sephadex™ G-25 spin columns 
(GE Healthcare, Germany) to remove free [y- 32 P]ATP. 
For probing PMM0685 mRNA or PMM0685, asRNA 
oligonucleotides PMM0685 rev and PMM0685-5'- 
end mRNA, respectively, were used (Supplementary 
Table SI). 

Transcript probes. PCR fragments containing a T7 recog- 
nition site were synthesized for in vitro transcription from 
the T7 promoter in the presence of [oc- 32 P]UTP (Hartmann 
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Analytic, Germany) using T7 MAXIscript (Ambion, 
USA). Sense and antisense probes were prepared for the 
following genes: PMM0684, PMM0685, PMED4 07411, 
PMED4_07431. Oligonucleotides used for template PCR 
amplification are listed in Supplementary Table SI. 

Northern blot analysis 

RNA samples (1.5 ug total RNA or 1.3 or 2.6pmol in vitro 
transcripts) were supplemented with loading buffer, 
denatured at 65°C for 5min, and then briefly chilled on 
ice before electrophoretic separation on 7M urea-6% 
polyacrylamide gels for 2h at 6 mA or 1.5% formalde- 
hyde-agarose or native gels for 2h at 80 V, respectively. 
Polyacrylamide gels were transferred to Hybond-N 
nylon membranes (Amersham, Germany) by electro- 
blotting for 1 h at 400 mA. Agarose gels were transferred 
overnight by capillary blotting using 20 x SSC (3 M NaCl, 
0.3 M sodium citrate, pH 7.0) as transfer solution. 
Membranes were hybridized overnight at 42 or 57° C 
with labelled oligonucleotide probes or labelled transcript 
probes, respectively, in 50% deionized formamide, 7% 
SDS, 250 mM NaCl and 120mM Na(P0 4 ) pH 7.2. 
Membranes were washed in 2x SSC/1% SDS for 
lOmin, in 1 x SSC/0.5% SDS for another lOmin and, ac- 
cording to signal intensity, 1-10 min in 0.1 xSSC/0.1% 
SDS. All wash steps were carried out at 5°C below the 
hybridization temperature with pre- warmed buffers. 
Membranes were exposed to imaging plates, and the re- 
sulting signals were visualized using a Personal Molecular 
Imager FX system with Quantity One software (Bio-Rad, 
Germany). 

Reverse transcription-PCR 

DNA was removed from total RNA using TURBO™ 
DNase (Ambion, USA). Treatment of RNA samples 
(3-1 5 ug) was performed in three consecutive incubation 
steps for 20 min each at 37°C using 6 U of DNase in total. 
The reaction was stopped with inactivation reagent 
(Ambion, USA), and the supernatant was purified and 
enriched by gel chromatography using RNA Clean & 
Concentrator M -25 according to the manufacturer's 
instructions (Zymo Research, USA). 

Reverse transcription (RT) was performed with 
0.75-1.5 jig total RNA using Superscript III (Invitrogen, 
Germany). Briefly, after denaturation at 95° C for 
3-5 min, a mixture containing RNA and 0.5 uM of 
target- and strand-specific oligonucleotide 'RT-mRNA' 
(for amplification of mRNA) or 'RT-asRNA' (for ampli- 
fication of asRNA) was incubated at 42°C for 20 min to 
permit annealing. The mixture was complemented with 
4U RiboLock RNase inhibitor (Fermentas, Germany), 
ImM dNTPs, lOmM DTT, 5mM MgCl 2 , 1 x reaction 
buffer (Invitrogen, Germany) and 200 U reverse transcript- 
ase in a total volume of 20 jil. For first-strand synthesis, the 
reaction was incubated at 42° C for 2h and the enzyme 
was inactivated for 5 min at 95°C. cDNA was amplified 
by PCR in a 50-|il reaction mixture comprising 0.2 uM 
oligonucleotides, 1 x HF reaction buffer (Finnzymes, 
USA), 1 U Phusion polymerase (Finnzymes, USA) and 
0.2 mM dNTPs. The following thermo profile was 



applied: 98°C for 30 s; 32 cycles of 98°C for 10 s, 54°C for 
30 s, and 72°C for 3 min; and 72°C for 5 min. To amplify 
cDNA derived from mRNA, the primers RT-mRNA 
(reverse primer) and GI-II-mRNA (forward primer) were 
used. To amplify cDNA derived from asRNA, the primers 
RT-asRNA (reverse primer) and RT-mRNA (forward 
primer) were used. Oligonucleotides used for RT-PCR 
are listed in Supplementary Table S 1 . The lack of contam- 
ination with genomic DNA was verified by omitting the 
enzyme in an additional reverse transcription control 
reaction. 



RESULTS 

Purification of Prochlorococcus RNase E 

Expression of recombinant Prochlorococcus MED4 
RNase E (PMM1501) protein required the adjustment 
of the codon usage of the first 63 nt and the last 42 nt of 
PMM1501 for sufficient expression in E. coli. C-terminally 
His-tagged pQE70-PMM1501 was overexpressed in E. coli 
strain Arctic Express (DE3) RIL and purified under 
native conditions by affinity column chromatography. 
Recombinant PMM1501 protein migrated according to 
its theoretical molecular weight of 70 kDa in SDS-poly- 
acrylamide gels (Figure 1 and Supplementary Figure SI). 
Mass spectrometry analysis of the purified protein fraction 
did not reveal any contamination with E. C6>//- specific 
RNase E protein or other host ribonucleases (data not 
shown). Furthermore, western blots of the purified 
protein fraction were immunodecorated with antibodies 
targeting E. C6>//- specific RNase E, RNase R or PNPase. 
None of these antibodies detected measurable amounts of 
the respective target protein in the elution fraction used 
for cleavage assays (Supplementary Figure S1A). To test 
for the presence of other sources of E. C6>//- specific exo- 
nuclease activity in the purified protein fraction, we 
in vitro transcribed a poly-A RNA, a substrate for exo- 
nucleases in general (RNase R, PNPase, RNase II) but not 
for RNase E, and incubated it with RNase R or recom- 
binant Prochlorococcus RNase E. Cleavage activity was 
only detected with RNase R but not with recombinant 
PMM1501 (Supplementary Figure SIB) verifying that 
the recombinant Prochlorococcus RNase E protein was 
not contaminated with nucleases from E. coli. 

To assess the quaternary structure of affinity-purified 
pQE70-PMM1501 protein, we analysed the oligomeriza- 
tion state by gel filtration. According to the calibration 
profile, the main part of purified recombinant MED4 
RNase E was eluted probably as a tetramer (Figure 1), 
which would suggest a similar conformer organ- 
ization as that present in E. coli and Synechocystis 
PCC6803. These results are in good agreement with the 
presence of the highly conserved Zn-link motif in the 
Prochlorococcus RNase E sequence. Callaghan et al. (31) 
demonstrated that the Zn-link is involved in the organ- 
ization of a functional active quaternary structure of 
RNase E/G in E. coli by mediating its dimerization and 
tetramerization. 
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Figure 1. Gel filtration profile of affinity-purified pQE70-PMM1501 protein. Sizes of RNase E conformers were estimated using Bio-Rad's gel 
filtration standard (black arrows) containing thyroglobulin (670 kDa), y-globulin (158 kDa), ovalbumin (44 kDa), myoglobin (17 kDa) and vitamin 
B12 (1.35 kDa). Elution peaks of RNase E conformers are indicated by black arrows. Peaks at elution fraction 15 and 21 correspond to buffer 
components. The inset on the right shows a 10% SDS-PAGE gel that was loaded with elution fractions of the gel filtration. 



Determination of the physicochemical properties of 
Prochlorococcus RNase E 

To establish optimal RNA cleavage conditions and to de- 
termine the physicochemical properties of MED4 RNase 
E, we assessed cleavage of a fTuorescently labelled synthet- 
ic oligonucleotide substrate containing the RNase E rec- 
ognition site of Synechocystis PCC6803 psbA leader (32) 
(see 'Material and Methods' section). We found that 
monovalent cations are not required whereas divalent 
cations are essential for MED4 RNase E enzyme activity 
with an optimum at 12mM Mg 2+ (Supplementary 
Figure S2). The enzyme was active between pH 5 and 
11, with optimal activity at pH 9 (Figure S2D). The 
optimal temperature for enzyme activity was 35°C 
(Supplementary Figure S2E), 10° above the optimal 
growth temperature for Prochlorococcus MED4 (33). 
Interestingly, the temperature optima for RNase E 
activity from Synechocystis PCC6803 [courtesy of M. 
Asayama (32)] was at 25°C (data not shown), which 
grows optimally at 30°C (34,35). Similar growth-deviant 
enzyme optima have been reported for purified glutamine 
synthetases of Prochlorococcus PCC9511 (a clonal culture 
of MED4) and Synechocystis PCC6803 with maximal 
enzyme activities at 55 and 35°C, respectively (36).The 
MED4 RNase E enzyme was quite stable up to 52° C, 
whereas a nearly complete inhibition of the enzyme was 
observed at 56°C (Supplementary Figure S2F). 

Transcript organization and expression of phage 
infection-induced transcripts located in genomic 
island II and in the ribosomal gene cluster 

A previous study on global expression profiles of 
Prochlorococcus MED4 transcripts during phage infection 



revealed the induction of a 1.5-kb region in genomic island 
II consisting of three small ORFs, from PMM0684 to 
PMM0686. Only PMM0686 has a proposed function 
(clpS-like protease) (27). Here, we investigated the 
genetic organization of this region in more detail. RT- 
PCRs and northern blot analysis revealed that this 
region is organized as a polycistron extending into two 
upstream genes, ssrS (6S RNA, Yfr7) and PMM0683 
(purK), and one downstream gene, yfr3 (Figure 2A and 
B). In total, the polycistronic region encompasses 3kb 
from nucleotide position 651 956-654 987, as determined 
from RT-PCR data and verified by northern blot analysis 
(Figure 2B and C). In another global microarray analysis, 
asRNAs partially covering the complementary strand of 
the same polycistronic region were identified (37). 
According to RT-PCR data and northern blot analysis, 
these asRNAs originate from one large precursor of 3.5 kb 
at nucleotide 651 469-654 987 on the complementary 
strand (Figure 2B and C), leading to complete coverage 
of the polycistronic region in genomic island II. Both 
RNAs on the sense and antisense strand are induced or 
stabilized during phage infection (Figure 2B and 
Supplementary Figure S3) (27,37) supporting our hypoth- 
esis that mRNA-asRNA duplex formation can protect 
RNA from RNase E degradation. Similar cleavage 
patterns were observed during exponential growth and 
phage infection despite changes in relative accumulation 
of the different bands (Figure 2B). Another region that is 
induced during phage infection is that of the ribosomal 
protein gene cluster. An increase of gene expression has 
been measured for rpsll (PMM1536), rpl6 (PMM1544) 
through rpl5 (PMM1546) and rpll6 (PMM1551) (27). 
These five genes belong to two adjacent operons ranging 
from PMM1533 to PMM1558 (38). Using specific primers 
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Figure 2. Organization of phage-induced transcripts from genomic island II. (A) Genome arrangement and transcript organization of the section of 
genomic island II starting with yfrl (46). Protein-coding genes and non-coding RNA genes (yfr) are shown in dark-blue and light-brown boxes, 
respectively. Transcript lengths as determined by RT-PCR are marked with horizontal arrows, blue for the forward strand and red for the reverse 
strand. The vertical arrows correspond to the PCR primers. Black arrows indicate primary transcript starts that were determined by 454 deep 
sequencing (unpublished data). Binding sites of probes used for northern hybridization shown in (B) are indicated by boxes and labelled with letter A 
in blue for the forward strand or B in red for the reverse strand. (B) Expression changes of PMED4_07411 in genomic island II and its asRNA 
during phage infection as determined from northern blot analysis. Lane C — RNA from uninfected control cells. Lanes 1 and 4 — RNA from cells 
during phage infection 1 or 4h after phage addition. Total RNA (1.5 jig) was loaded per lane. A blot with higher resolution of a section of 
asPMED4_07411 is presented to the right of the asRNA hybridization. As a loading control rnpB was hybridized. Transcript lengths were estimated 
from Fermentas RNA high-range marker (M). Probes used for northern hybridizations are indicated by boxes labelled with letter A in blue for the 
forward strand or B in red for the reverse strand. (C) RT-PCRs of genomic island II (GI II) region on the forward (fw) and reverse (rv) strand 
revealing an asRNA of ~3.5kb. (D) RT-PCR of the asRNA transcript of ribosomal protein cluster (asrps/1) revealing an asRNA of approximately 
7kb. PCR fragment lengths were estimated from lambda PstI marker (M). The identity of cDNAs was verified by Sanger sequencing of PCR 
products. (E) Northern blot hybridized with a probe targeting the antisense region of rpsll. The asrpsll probe was synthesized using oligonucleotides 
asrpsll_FW and asrpsll_RV (Supplementary Table SI). M — fermentas RNA high-range marker. 



in RT-PCR, we detected a 7-kb asRNA from nucleotide 
1 476 365-1 483 393 that covers the complete 
phage-induced region (Figure 2D and E). Notably, the 
detection of a 7-kb asRNA is the first report of a prokary- 
otic asRNA with such an extraordinary length. However, 
asRNAs covering more than one gene have been reported 
before in the bacterium Listeria monocytogenes (39). 

RNA-RNA duplex formation masks recognition sites 
of RNase E 

We established an in vitro RNase E cleavage assay to 
follow up our initial questions about how the activation 



of genes in genomic island II is regulated during phage 
infection and how these RNAs are protected from degrad- 
ation in light of the enhanced transcription of RNase E. 
The mRNAs of four small ORFs of genomic island II 
(PMM0684, PMM0685, PMED4 07411, PMED4_7431), 
which were induced during phage infection (Figure 2B 
and Supplementary Figure S3), were synthesized in vitro 
and incubated with purified recombinant MED4 RNase 
E. All four transcripts were recognized by RNase E 
and were susceptible to specific cleavage (Figure 3). 
RNase E and RNase III cleavage sites of PMM0685 
were investigated in more detail by northern hybridization 



4896 Nucleic Acids Research, 2011, Vol 39, No. 11 



PMM0684 mRNA PMED4_07411 mRNA PMM0685 mRNA PMED4_07431 mRNA 




Figure 3. RNase E cleavage assay of four in vitro transcripts of genomic island II with (+) and without (— ) recombinant Prochlorococcus RNase E. 
Cleavage fragments were separated on a 7M urea-6% PAA gel and stained with ethidium bromide. Fragment sizes were estimated from an NEB 
ssRNA marker. 



(Supplementary Figure S4) using specific oligonucleotide 
probes targeting the 5'- and Spends, respectively, of the 
ribonuclease-treated transcripts. The RNase E cleavage 
sites identified are listed in Table 1. All four cleavage 
sites were located in loop regions, indicating that the sec- 
ondary structure is more relevant than sequence specificity 
for RNase E recognition (Figure 4B). RNase III cleavage 
sites were located in double-stranded RNA regions 
(Figure 4B) and are specified in Table 2. Figure 4 A 
shows the results of RNase E, RNase III and RNase Tl 
digestions for single-stranded and duplex RNAs of 
PMM0685 in v/Yrosynthesized RNA. Both sense and anti- 
sense RNAs of duplex RNAs had identical sizes and an 
exact overlap. Incubation of single-stranded PMM0685 
RNA with RNase E or RNase III yielded specific 
cleavage products, whereas RNase Tl, which cleaves 
single-stranded RNA after guanine residues, led to 
complete digestion of PMM0685 in vitro- synthesized 
RNA. In contrast, the a priori duplex formation of sense 
and antisense RNA before RNase E treatment prevented 
RNase E digestion. RNase III treatment of duplex RNA 
resulted in complete digestion of the duplex, whereas 
treatment with RNase Tl did not show any cleavage 
products. These findings suggest that the secondary and 
tertiary structures of both RNA molecules were altered 
and hybridized to each other over their full-length. 
Clearly, the resultant structural changes hinder substrate 
accessibility of ribonucleases that recognize single- 
stranded regions (RNase E and RNase Tl) and promote 
enzyme activity of RNA double strand-dependent 
enzymes like RNase III. 

RNA duplexes that overlapped only partially with over- 
hangs at the 5'- or 3'-end lead to partial degradation of 
the duplex by RNase E (Figure 5). As expected, digestion 
patterns indicated that double-stranded regions were 
still protected, whereas single-stranded regions became 
susceptible to RNase E cleavage (Figure 5). When a 



Table 1. RNase E cleavage sites of PMM0685 in vitro transcript 



Position within RNA backbone (5'-3') 


Cleavage site (5-3') 


53, 54 


U GAAAA 


80, 81 


AU UACU 


153, 154 


AAG CUU 


162, 163 


CU CAAG 



3-kb polysictronic mRNA in vitro transcript was 
annealed with a 3.5 kb in vitro asRNA (that mimic the 
in vivo transcripts) complete protection of the polycistron 
after RNase E treatment was observed (Supplementary 
Figure S5). The in vitro assay suggests that indeed the 
3.5 kb asRNA located in genomic island II is important 
for the stabilization of mRNAs located opposite the 
asRNA. These results suggest that complete coverage of 
the mRNA by an asRNA is required to prevent RNase E 
digestion. It is conceivable that the same protective mech- 
anism applies for the ribosomal protein gene cluster as 
the asRNA covers the complete phage induced region 
(Figure 2). Therefore, complete overlap of duplex RNAs 
may be a general requirement for mRNA protection. 



DISCUSSION 

In this study, we describe a new function for asRNA, a 
class of regulatory RNA that has gained considerable at- 
tention due to its recently discovered ubiquitous occur- 
rence in all kingdoms of life. Previously, we (27) 
suggested that the subset of host genes upregulated 
during phage infection might be regulated by the phage 
or might be part of a host-specific response. For either 
scenario, it remained unclear how asRNAs are regulated. 
A prolongation of mRNA half-life initially appeared 
unlikely because RNase E expression was stimulated as 
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ss RNA 



duplex RNA 



Table 2. RNase III cleavage sites of PMM0685 in vitro transcript 



nt 

1000 

500 

300 



150 



80 



50 



M Ctrl oc 



z z 



Ctrl 



z 
oc 



z 
oc 



ss RNA 



— duplex RNA 




I: internal loop 
E: external loop 
S: stem 



Figure 4. Ribonuclease cleavage assay of 1.3pmol single-stranded and 
2.6pmol duplex PMM0685 in vitro RNA. (A) Ethidium bromide- 
stained 7 M urea-6%PAA gel of PMM0685 single-stranded (ss) and 
duplex in W/ro-synthesized RNA (ctrl) treated with RNase E, RNase III 
or RNase Tl. RNA was incubated with RNases for 15min. M denotes 
ssRNA marker from NEB. Though gel electrophoresis was performed 
under stringent denaturing conditions, only small amounts of duplex 
RNA were denatured and migrated according to the single-stranded 
RNA population. (B) The secondary-structure of PMM0685 was 
folded with RNAfold and drawn using VARNA. RNase E and 
RNase III cleavage sites, which were deduced from northern blot hy- 
bridizations with two oligonucleotide probes targeting the 5'- or 3'-end 
of PMM0685 in vz/rosynthesized RNA, are marked with red arrows or 
blue arrows, respectively. 



well. Here, we demonstrate the asRNA-mediated regula- 
tion of RNA stability, directly influencing the half-life of 
mRNAs that are protected from RNA decay by asRNAs 
that completely mask RNase E recognition sites. Whether 
the increase in transcript level is a final defensive step of 
the host during phage infection or whether increased 
mRNA levels are supportive of viral replication [in a 
similar fashion as has been suggested for phage-encoded 
Dl protein (40)] remains to be investigated in future 
studies. The asRNA-modulated regulation of gene expres- 
sion during phage infection might be a common feature in 



Position within RNA 
backbone (5 f -3 f ) 



Cleavage site (5'-3 r ) 



16,17/27,28 




C C/G G 


47,48/225,226 




C A/G U 


68,69/220,221 




C A/G U 


72,73/216,217 




G_A/C_U 


77,78/86,87 




G C/C G 


89,90/201,202 




C U/G A 


97,98/193,194 




C G/G U 


102,103/187,188 




C A/G U 


110,111/136,137 




U A/G U 


115,116/131,132 




C A/G U 


141,142/183,184 




C U/G G 


160,161/166,167 




C_C/G_G 


A 


ss RNA 


duplex RNA 
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> 
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Figure 5. Ribonuclease cleavage assay of 1.3pmol single-stranded and 
2.6pmol partially overlapping duplex PMM0684 in vzYro-synthesized 
RNAs. (A) 7M urea-6% PAA gel of PMM0684 single-stranded and 
duplex in vitro RNAs (ctrl) treated with RNase E, RNase III or RNase 
Tl. M denotes ssRNA marker from NEB. (B) Scheme of duplex for- 
mation of PMM0684 in W/ro-synthesized mRNA (blue) and asRNA 
(red). Regions of duplex RNA remaining after digestion with RNase 
E or RNase III in (A) are labelled with corresponding numbers in (B). 



both the host and the phage. Millard et al (41) reported 
an asRNA that physically connects between the photosyn- 
thesis gene psbA encoding the Dl core protein of photo- 
system II with a homing endonuclease (F-Cphl) in the 
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cyanophage S-PM2. Transcripts of the asRNA as well as 
of psbA and F-Cphl are induced during the lytic infection 
cycle of S-PM2 in its host Syne cho coccus sp. WH7803. 

RNase E activity is known to impact the stability of 
both host and phage mRNA during infection of E. coli 
by bacteriophages T7 and T4 (42,43). In both cases, the 
C-terminal structural domain of RNase E is involved. For 
example, the T7 protein kinase phosphorylates the struc- 
tural domain of host RNase E. However, this region is not 
present in the Prochloro coccus MED4 homologue. 
Therefore the indirect regulation of RNase E activity in 
Prochlorococcus by RNA duplex formation reported here 
reflects an alternative mechanism of regulation of enzyme 
activity that appears to be required for an organism 
lacking the region modified in other host-phage systems. 

Our results clearly show that the duplex RNAs resulting 
from mRNA-asRNA hybridization are perfect substrates 
for RNase III in vitro (Figure 4). Therefore, an interesting 
question is how these duplex RNAs become protected 
from RNase III degradation in vivo. Nothing is known 
about RNase III protein stability in Prochlorococcus; 
however, the RNase III transcript, along with the 
majority of the transcriptome, is down-regulated during 
phage infection (27). Furthermore, ribosome activity on 
the mRNAs is likely to interfere with the full-length hy- 
bridization of both transcripts in vivo. In addition, other 
factors might exist in the cell that protect long duplex 
complexes from degradation. 

The modulation of RNase E activity by RNA-RNA 
interaction has been reported before. However, base 
pairing of the small regulatory RNA RyhB with the 
S'-UTR of sodB mRNA, changing the structure of the 
mRNA leader, results in an additional RNase E 
cleavage site, which in turn triggers another pathway for 
sodB RNA decay (44). Hence, in this example, RNA- 
RNA interaction does not hide but produces new RNase 
E recognition sites. The RNase E-dependent regulation of 
sodB turnover is very complex, and enzyme activity is also 
influenced by RNA-protein interaction between sodB and 
Hfq (an RNA chaperone that recognizes single-stranded 
AU-rich regions in the same manner as RNase E). Both 
RNase E and Hfq can bind mRNAs close to each other, 
resulting in protection of the RNase E substrate from 
nucleolytic degradation (44). Such a complex regulation 
cannot be expected for Prochlorococcus MED4 because an 
hfq homologue is lacking in the genome. Of the 12 com- 
pletely sequenced Prochlorococcus strains, only MIT9313 
and MIT9303, both belonging to the evolutionarily oldest 
Prochlorococcus clade, possess an hfq homologue. 

Alternatively, RNase E activity can also be influenced 
by inhibitory RNase-binding proteins that dynamically 
regulate RNA decay and processing. In Syne cho coccus 
elangatus PCC 7942, the two heat shock proteins Dnak2 
and DnaJ2 inhibit RNase E activity in an ATP-dependent 
manner, suggesting that both proteins are involved in 
RNA degradation through interaction with RNase E 
(45). Both proteins are present in Prochlorococcus and 
hence could provide an additional regulatory circuit for 
RNA half-lives. 

The reduced suite of protein regulators and the dispro- 
portionately high number of regulatory RNAs in 



Prochlorococcus has been linked to its adaptation to an 
extremely oligotrophic habitat, combined with a 
streamlined genome (37). This study provides the first 
insights into the functional role of asRNA regulators 
and emphasizes the importance of this class of regulatory 
molecules in the global network of gene regulation. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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